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1. INTRODUCTION 

Clarification of the nature of the finite temper- 
ature QCD matter and its temperature driven 
phase transition or crossover is one of the ma- 
jor goals of QCD calculations on the lattice. For 
simplicity, I use the term "transition" for both 
genuine phase transitions and sharp crossovers in 
this report, unless explicitly specified. In a study 
of the order of a transition, the scaling analysis is 
a powerful tool. Unlike for the SU(3) pure gauge 
theory where the finite spatial size of the lattice is 
the only cause to smoothen the singularity at the 
transition point [1], near the chiral transition in 
QCD the quark mass is expected to play a dom- 
inant role provided that the spatial lattice size is 
not too small [2]. The scaling analysis is particu- 
larly important in a study of the chiral transition 
with two flavors because theoretical studies do 
not have a definite prediction about the order of 
the transition in this case [3,4]. This year, new 
studies along this line were presented. 

Most numerical studies in finite temperature 
QCD have been performed with staggered quarks. 
Although the flavor structure is much simpler 
with Wilson quarks, the violation of chiral sym- 
metry due to the Wilson term poses several con- 
ceptual and technical difficulties. Recent efforts 
with Wilson quarks revealed several unexpected 
properties and a more complicated phase struc- 
ture than that for staggered quarks [5]. In order 
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to make a reliable prediction for the continuum 
limit, however, it is essential to remove the de- 
pendence on the lattice formulation of quarks in 
the results. Several fundamental progresses were 
reported for Wilson quarks at the conference. 

In a more realistic case including the strange 
quark, the status is still not clear enough. With 
staggered quarks, the Columbia group found a 
time history that suggests a crossover at a point 
where m a is smaller than its physical value and 
Tii u and are larger than their physical values 
[6]. Because the results for two flavors suggest a 
crossover at m a = oo and at these values of m u 
and rrid, this implies that the transition in the real 
world is also a crossover provided that the depen- 
dence on the light quark mass is negligible. On 
the other hand, a recent study of an effective o 
model suggests that the nature of the transition 
can sensitively depend on the light quark mass 
[4,7]. Therefore, it is important to determine the 
global phase structure in the full coupling param- 
eter space. A new result with Wilson quarks was 
presented. 

For the nature of the finite temperature QCD 
matter, simulations on larger lattices were made 
to determine the transition temperature and the 
equation of state; the behavior of high temper- 
ature gluon propagators was studied; quenched 
calculations on the finite temperature hadronic 
matter were performed; and properties of a va- 
lence quark chiral condensate were studied. 

QCD at finite density is another highly impor- 
tant subject especially in phenomenological appli- 
cations to heavy ion collisions and astrophysical 
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Table 1 

Critical exponents of the three dimensional 0(4) [10] and 0(2) [11] Heisenberg models, the meanfield 



theory, and Np — 


2 QCD with stag 


gered quarks [12]. 










0(4) 


0(2) 


MF 


Karsch-Laermann 


1/(36 


0.537(7) 


0.602(2) 


2/3 


0.77(14) 


ft 


1/6 


0.2061(9) 


0.2072(3) 


1/3 


0.21 - 0.26 


chiral cumulant 


1 - 1/6 


0.7939(9) 


0.7928(3) 


2/3 


0.79(4) 


Xm 


{1-0)1158 


0.331(7) 


0.394(2) 


1/3 


0.65(7) 


Xt 


a/156 


-0.13(3) 


-0.003(4) 





-0.07 - +0.34 


specific heat 



processes. Nevertheless, a reliable numerical sim- 
ulation is still difficult as a consequence of the 
complex action problem [8]. A few new studies 
were reported this year. 

In this report I attempt to review these new 
developments. In Sect. 2, I review recent studies 
on the chiral transition with two flavors of stag- 
gered and Wilson quarks. Sect. 3 deals with the 
influence of the strange quark. In Sect. 4, we 
concentrate on other subjects in the QCD ther- 
modynamics. Recent progress at finite chemical 
potential is discussed in Sect. 5. A brief sum- 
mary is given in Sect. 6. For the previous status 
of investigations, I refer the reader to Refs. [9,5,7]. 

2. CHIRAL TRANSITION IN QCD 
WITH TWO FLAVORS 

Understanding the nature of the finite tem- 
perature transition in QCD with two degenerate 
light quarks {Np = 2) is an important step to- 
ward the clarification of the transition in the real 
world. Based on a study of an effective a model 
Pisarski and Wilczek discussed that the transi- 
tion in the chiral limit (chiral transition) of QCD 
with two flavors is either of second order or of 
first order depending on the strength of the Ua(1) 
anomaly term at the transition temperature [3]. 
Three years ago, Wilczek and Rajagopal argued 
that, in case that the chiral transition is of sec- 
ond order, QCD with two flavors will belong to 
the same universality class as the three dimen- 
sional 0(4) Heisenberg model [4]. This provides 
us with several useful scaling properties that al- 
low us to make clear the order of the transition. 
Because no known lattice fermions have the full 
chiral symmetry on finite lattices, the appearance 
of the 0(4) scaling will be a useful touchstone to 



test the recovery of the full chiral symmetry on 
the lattice if the chiral transition is of second or- 
der. 

Scaling properties are described in terms of 
critical exponents. Consider a spin model at tem- 
perature T near the transition temperature T c 
and at small external magnetic field h. The ex- 
ponents a, (5, 7, and 8 are defined for the specific 
heat C , the magnetization M, and the magnetic 
susceptibility x by 



C(t,h = 0) 
M{t,h = 0) 

X(t,h = 0) 
M{t = 0,h) 



\t\~ a 

\tf (t<0) 
t-< (t > 0) 



(1) 

(2) 
(3) 
(4) 



where t = [T -T c {h = 0)]/T c {h = 0) is the reduced 
temperature. These critical exponents satisfy the 
following two scaling relations: a = 2 — (5{8 + 1) 
and 7 = (5{8 — 1), so that only two of the four ex- 
ponents are independent. Some 0(4) exponents 
are listed in Table 1. In QCD, h corresponds to 
the quark mass and M corresponds to the chiral 
condensate. 

2.1. Scaling study with staggered quarks 

The conjecture by Wilczek and Rajagopal has 
been tested first for staggered quarks [12]. Nu- 
merical results obtained so far for the transition 
with Np = 2 staggered quarks are consistent with 
a second order chiral transition: No signs of dis- 
continuities are found at m q a > 0.0125 for 7V t = 4 
and 6, where 7V t is the lattice extension in the 
temporal direction [6,13,14]. Also a finite size 
scaling for a first order transition is not compati- 
ble with the data obtained at m q a = 0.01 - 0.0125 
and N t = 4 [14,15]. 

The last two columns in Table 1 show the nu- 
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Figure 1. Pseudocritical coupling f3 c versus m q /T 
at N t = 4 [16]. This figure is an update of 
Fig. 6 in Ref.[12] including new results obtained 
on 12 3 x 4 and 16 3 x 4 lattices (a part of diamonds 
at m q /T = 0.08, 0.13, 0.3 and 0.6). 

merical results for critical exponents by Karsch 
and Laermann obtained at m q a = 0.02 - 0.075 
on an 8 3 X 4 lattice [12]. The result for 6 is con- 
sistent with the 0(4) value, while the results for 
a and f3 are in disagreement with the 0(4) ex- 
ponents. Clearly, we need data with better accu- 
racy (on larger lattices and at smaller m q a). The 
Bielefeld group continued their efforts to improve 
the simulation and preliminary results on larger 
spatial lattices (12 3 and 16 3 ) were presented at 
the conference [16]. The results obtained so far 
are consistent with those on the 8 3 x 4 lattice (cf. 
Fig. 1). Another test of the 0(4) scaling per- 
formed by DeTar [9] will be discussed later. 

With staggered quarks, however, several 
caveats are in order because Np = 2 staggered 
quarks are realized on the lattice by introducing 
a fractional exponent to the fermionic determi- 
nant of Np = 4 staggered quarks. Therefore, 
(i) the symmetry in the chiral limit on a lattice 
with finite lattice spacing is the 0(2) symmetry 
of Np = 4 staggered quarks [17], and not the 
0(4) symmetry we expect to find in the contin- 
uum limit, and (ii) the action is not local. The 
correct continuum chiral limit with the 0(4) sym- 
metry will be obtained only when we first take the 



continuum limit (3 — > oo and then take the chiral 
limit. Choosing a too small m q compared with 
the lattice spacing may lead us either to wrong 
0(2) exponents, or to some non-universal behav- 
ior due to the lack of locality. 

Another theoretical possibility (both for stag- 
gered and Wilson quarks) is the appearance of 
the meanfield (MF) exponents [18]. Out of the 
critical region (Ginzburg region), the order pa- 
rameter may be determined by a MF theory [7]. 
For example, the critical region in the BCS theory 
is known to be quite narrow so that many prop- 
erties of superconductors can be calculated with 
the MF theory by Ginzburg and Landau. In 3D 
and 4D Gross-Neveu models, Kocic and Kogut 
reported the absence of critical regions [18]. Nar- 
rowness in the width of the critical region in the 
BCS theory is closely connected to the fact that 
Cooper pairs are very loosely bounded. In QCD, 
because quarks are tightly bounded at low tem- 
peratures, we may expect the critical region to be 
less narrow. 

These possible exponents are listed in Table 1. 
It is important that we are able to make a nu- 
merical distinction between these different sets of 
critical exponents. Unfortunately, 0(4) and 0(2) 
exponents are very close to each other. 

2.2. Phase structure and problems with 
Wilson quarks 

In the case of Wilson quarks, the lattice action 
lacks chiral symmetry due to the Wilson term. 
This introduces possible O(a) effects to physical 
quantities and relations characteristic for chiral 
symmetry. However, when we are close enough 
to the continuum limit, we expect that the sym- 
metry breaking effects become sufficiently weak 
so that we see the 0(4) scaling when the chiral 
transition is of second order. In this case the ap- 
pearance of the 0(4) scaling is a useful test to see 
if the effects from the Wilson term are small. 

Although chiral symmetry is broken explic- 
itly, numerical simulations show that the relation 
m\ oc m q at small m q is well satisfied in the low 
temperature phase for all values of (3, where m q 
is defined via an axial Ward identity [19-21]. It 
is therefore natural to define the chiral limit K c 
by m T = or m q = at T = 0. The loca- 
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Figure 2. Phase diagram for Np = 2 QCD with 
Wilson quarks using the standard action. See text 
for details. Lines are to guide the eyes. 

tion of K c for Np = 2 is shown by the dark-gray 
line in Fig. 2 collecting recent data from various 
groups [5]. There exist small discrepancies (at 
most of the order of 0.01) between the values for 
K c from m x = and m q = which are attributed 
to O(a) effects due to the Wilson term, and also 
to the errors from the extrapolation of and 
m q in 1/K. The light-gray lines K t (N t ) indicate 
the finite temperature transition determined from 
rapid changes in physical observables. Because 
the line K c is the renormalization group flow line 
for massless QCD, the finite temperature transi- 
tion for massless QCD is located at the crossing 
point of K t and K c . We shall therefore denote 
this point as the chiral transition point /3 ct . No- 
tice that there are possible small ambiguities of 
0(a) in the location of the K c line and the point 

Art- 

In the high temperature phase, we encounter 
no singularities even just on the K c line, where 
m q vanishes [21] but the finite temperature pion 
screening mass m T (T > 0) remains non-zero 
[22,23]. Also with staggered quarks, we expect 
a similar absence of singularity on the m q = 
line in the high temperature phase. In the low 
temperature phase, on the other hand, there ap- 
pear small eigenvalues of the quark matrix near 
K c which make the number of CG iterations re- 



quired to invert the quark matrix to a given ac- 
curacy, very large [24]. The Tsukuba group per- 
formed simulations on the K c line by reducing (3 
from the high temperature side, and looked for a 
point where singularity first appears in the sense 
that the number of CG iterations becomes very 
large. They found that the location of the point 
is given by (3 ~ 3.9 - 4.0 for N t = 4 and 4.0 - 4.2 
for Nt = 6. These two points are consistent with 
the ones obtained from linear extrapolations of 
the K t lines [23]. Therefore, they identified these 
points with f3 c t(N t ). 

Recently, Aoki, Ukawa and Umemura stressed 
that, in understanding the phase structure of Wil- 
son quarks at finite temperatures, it is important 
to clarify the meaning of the critical line at T > 
[25]. They defined the finite temperature K c line 
(which I will denote in the following as the Kf >0 
line) by m^(T > 0) = 0, and studied it in the 
parameter space including the region above K c . 
Let us denote K^ >0 at T = as K^=° . Note 
that the K^=° line defined by m^(T = 0) = has 
no width and locates in the K c band shown in 
Fig. 2 whose width is cause by O(a) ambiguities 
mentioned above. 

Motivated by their analytic results of the fi- 
nite temperature 2D Gross-Neveu model. Aoki 
et al. performed simulations of Np = 2 QCD 
on an 8 3 x 4 lattice and found that the K^ >0 
line, which first flows along the K c line in the 
low temperature phase, sharply turns back up- 
ward in the vicinity of (3 = 4.0 (cf. the dashed 
line in Fig. 2). The study of the Gross-Neveu 
model suggests that the lower part of the K^ >0 
line locates slightly (probably 0(a) « 0(l/N t )) 
above the Kf =0 line. If this is the case, m T (T c > 
T>0) « 0(1/N t ) > on the K^=° line and will 
vanish only in the continuum limit. The study of 
the Gross-Neveu model also suggests that, when 
Nt is increased, the turning point of the Kf >0 
line moves along the K c line toward larger (3 and 
touch the weak coupling chiral point K = 1/8 in 
the limit Nt = oo. The region surrounded by the 
Kf >0 line is the parity broken phase proposed by 
Aoki for T = QCD [26]. Concerning the rela- 
tion between Kf >0 and K t in QCD, they argued 
that the K t line crosses the Kf =0 line and runs 
past the turning point of the Kf >0 line without 
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touching this line (cf. Fig. 4 in Ref.[25]). 

The shape of the K^ >0 line shown in Fig. 2 
is in accordance with the absence of singularity 
in the high temperature phase discussed above. 
The location of the turning point is consistent 
with /3 ct estimated by the Tsukuba group. The 
suggested O(a) gap between the K^ >0 and K^ =0 
lines in the low temperature phase was previously 
observed by the Tsukuba group [23]. We note 
that the possibility is not excluded that the K t 
line touches the turning point of Kf >0 . In this 
case, the crossover K t becomes a genuine transi- 
tion at that point. 

When we view the lattice QCD at finite 7V t and 
finite a as a statistical system, the phase structure 
is understood in terms of the critical line K^ >0 . 
On the other hand, because our final goal is to 
investigate the chiral transition of QCD in the 
continuum limit, it is important to know the rela- 
tion between the K c line and the K t line and also 
to trace the recovery of chiral symmetry which is 
broken on the lattice. This leads us to the K c line 
and (3 ct which have O(a) ambiguities as discussed 
above. In order to recover the finite temperature 
chiral transition in the continuum limit, the gaps 
between K^ >0 and ifj =0 and between K t and 
the turning point of K^ >0 should be also at most 
of O(a). The data actually show that these gaps 
are small to the present numerical accuracy. The 
continuum limit is not affected by these details. 

The important progress made by Aoki et al. is 
that the existence of the points where m^{T > 
0) = is explicitly shown in a toy model, and 
that the phase structure inferred from this model 
is confirmed in QCD at least partly. This provides 
us with a more rigid theoretical basis to extract 
the properties of the chiral transition in the con- 
tinuum limit from physical observables near /3 ct . 

From the behavior of physical quantities on the 
K c line, the chiral transition at (3 = (3 ct is found to 
be continuous both for N t = 4 and 6 as in the case 
of staggered quarks [23]. Nevertheless, it seems 
less hopeful to observe the desired 0(4) scaling on 
these lattices. The main reason is the fact that 
the dependence of the transition on (3 and the 
quark mass are completely different from those 
we expect when the chiral transition is of second 
order [22,27,28]: Near the continuum limit, we ex- 



pect that, as the quark mass increases from the 
chiral limit, the transition becomes weaker with 
the quark mass and it becomes strong again when 
the quark mass is heavy enough to recover the 
first order transition of the SU(3) gauge theory. 
Contrary to this expectation, the MILC collabo- 
ration found that, when we decrease K from K c 
on an 7V t = 4 lattice, the transition on K t be- 
comes once very strong at K ~ 0.18 and becomes 
weaker again at smaller K [22]. On a lattice with 
N t = 6 they even found a first order transition at 
K = 0.17 - 0.19 [27]. 

Looking at the phase diagram shown in Fig. 2 
more closely, we note that these points of strong 
transition are just the ones where the K t lines 
become once very close to the K c line due to 
the sharp bend of the K c line at (3 ~ 5.0, which 
is caused by the cross-over phenomenon between 
weak and strong coupling regions of QCD. There- 
fore, it seems plausible that the strong transition 
is a result of lattice artifacts caused by this un- 
usual relation of the K t and K c lines [5]. The 
quark mass m q also shows an unexpected 7V t 
dependence in the high temperature phase at 
/3 < 5.3 [22,28], which is due to lattice artifacts, 
too. 

A naive way out of these lattice artifacts is to 
increase 7V t so that we have the chiral transition 
in the weak coupling region, f3 ct > 5.5. However, 
a previous study on an 18 3 x 24 lattice suggests 
that this requires 7V t > 18, although the spatial 
lattice size is not large enough [23]. 

2.3. Wilson quarks with improved action 

Problems discussed in the last part of the previ- 
ous subsection with Wilson quarks using the stan- 
dard action have motivated the Tsukuba group to 
study finite temperature physics using improved 
actions [28,29]. They studied the renormaliza- 
tion group improved gauge action proposed by 
Iwasaki [30] and combined it with the standard 
Wilson quark action. The phase diagram is given 
in Fig. 3. Unlike for the standard action (cf. 
Fig. 2), the distance between K t and K c grows 
now monotonically when we increase (3 from the 
chiral transition point [3 ct ~ 1.4. Correspond- 
ingly, the K t transition becomes rapidly weak as 
we increase (3 starting from (3 ct (cf. Fig. 4). The 
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Figure 3. Phase diagram for Np = 2 QCD with 
Wilson quarks coupled to a RG improved gauge 
action [29]. K c is the chiral limit determined by 
ra\ = on an 8 4 lattice and K t is the finite tem- 
perature transition line obtained on an 8 3 x 4 lat- 
tice. Lines are to guide the eyes. 



straight line envelop of m\ for N t = 4 shown in 
Fig. 4b coincides with m\ for N t = 8 and cor- 
responds to the chiral behavior ra\ oc m q in the 
low temperature phase. Deviation from this line 
signals the crossover of the system to the high 
temperature phase. The smoothness of physical 
observables strongly suggests that the transition 
is a crossover at p > f3 ct . We can also see that m\ 
on the K c line monotonically decreases to zero as 
P ~^ Pet + [29], implying that the chiral tran- 
sition is continuous. It has also been shown that 
the lattice artifact in m q in the high temperature 
phase mentioned before disappears with this im- 
proved action. 

These nice properties which are in accordance 
with naive expectations encouraged the Tsukuba 
group to begin a scaling study with Wilson 
quarks. A key quantity to study scaling proper- 
ties is the chiral condensate (magnetization). The 
naive definition of (^^) is not so useful because 
the chiral symmetry is explicitly broken due to 
the presence of the Wilson term. A proper sub- 
traction and a renormalization are required to ob- 
tain the correct continuum limit. Just like m q 
[19,20], a properly subtracted (^^) can be de- 



Polyakov 
improved 




0.5 1 
1/K-l/K 



Figure 4. The Polyakov loop and the pion screen- 
ing mass with a RG improved action at various f3 
on an 8 3 X 4 lattice [29]. 



fined via an axial Ward identity [19]: 
<**) su b = 2m q aZY / {Tr(x)ir(0)) 



(5) 



where Z is the renormalization coefficient. This 
definition is consistent with the identification 
of the magnetization by Rajagopal and Wilczek 
[4]. ( 1 5 , ^ , )sub was shown to have a non-vanishing 
value in the chiral limit in the confining phase of 
quenched QCD [31] and QED [32]. It is easy to 
see that Eq.(5) holds also with improved gauge 
actions by replacing Z. For our purpose, it is 
enough to use the tree value: Z = (2K) 2 . 

We now study {^ , ^ , ) su b at finite temperature. 
Results for N t = 4 with the improved action are 
shown in Fig. 5. From the universality argument 
we expect that magnetization M near the tran- 
sition point can be described by a single scaling 
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0.5 



1.5 



Figure 5. Subtracted chiral condensate {^ , ^ , ) su b 
as a function of 2m q a for Wilson quarks with a 
RG improved action on an 8 3 x 4 lattice [29]. 



function: 

M = h}/ 6 f(t/h}/P 6 ) 



(6) 
and 



where h is the external magnetic field (oc m q 
t is the reduced temperature. DeTar tested this 
scaling for the case of two flavor staggered quarks 
[9] and found that data are consistent with 0(4) 
and 0(2) scaling. 

Fig. 6 shows the result for the scaling function 
f{t/h}^ 6 ) = (**) sub //i 1/ ' 5 with the identifica- 
tion M = (**) su b, h = 2m q a and t = (3 - (3 ct . 
With fixing the exponents to 0(4) and MF val- 
ues given in Table 1, we adjust (3 ct to obtain 
the best fit. With the 0(4) exponents, we find 
that the scaling ansatz works remarkably well as 
shown in Fig. 6a. The resulting f3 ct = 1.34(3) 
is slightly smaller than ~ 1.4 which is the value 
obtained by a linear extrapolation of the K t line 
(cf. Fig. 3) and of rn\ curves on K c and K t [29]. 
However, the 0(4) universality predicts [4] that 
(3c{m q ) ~ fid oc m\ /l3S and m\ ~ (fi c - (3 ct )' y on 
K c and K t with 7 ~ 1.4, i.e. these lines should 
bend slightly near the chiral transition point to 
give a smaller (3 ct . Therefore, we conclude that 
f3 ct = 1.34(3) is consistent with the data. On the 
other hand, fixing the exponents to the MF values 
makes the fit less beautiful. The best fit shown in 
Fig. 6b is obtained with f3 ct = 1.48(3). Contrary 
to the case of the 0(4) exponents, (3 ct = 1.48(3) 
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Figure 6. Best fits for the scaling function with 
0(4) and MF exponents. The plot contains all 
data of Fig. 5 within the range < 2m q a < 0.9 
and P < 2.0. 



is too large to be accepted. 

In conclusion, the data obtained with the im- 
proved action are consistent with the predicted 
0(4) scaling but inconsistent with the MF scal- 
ing. The success of this scaling test strongly sug- 
gests that the chiral transition is of second or- 
der in the continuum limit. It also indicates that, 
with the improved action, the chiral violation due 
to the Wilson term is sufficiently weaker than that 
introduced by the non- vanishing m q at least for 
m q s studied here. For more confirmation of these 
results, direct extraction of each critical exponent 
is required. 
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Figure 7. Map of expected nature of the QCD 
transition for Np =2 + 1 QCD as a function of 
the u and d quark mass m u d and the s quark mass 



3. INFLUENCE OF THE STRANGE 
QUARK 

For QCD with Np > 3, Pisarski and Wilczek 
predicted a first order chiral transition [3]. This 
was confirmed by numerical simulations both 
with staggered [6,33] and with Wilson quarks 
[23]. Off the chiral limit, the first order tran- 
sition smoothens into a crossover at sufficiently 
large m q . Therefore, the nature of the transition 
sensitively depends on Np and m q . This means 
that, in order to study the nature of the tran- 
sition in the real world, we should include the s 
quark properly whose mass m s is of the same or- 
der of magnitude as the transition temperature 
T c ~ 100 - 200 MeV. 

3.1. Expected phase structure 

Following Brown et al. [6], we summarize in 
Fig. 7 what we expect about the nature of the fi- 
nite temperature transition as a function of quark 
masses neglecting the mass difference among u 
and d quarks (Np = 2 + 1). When all quarks 
are heavy, the transition is of first order as ob- 
served in the SU(3) pure gauge theory. The limit 
m 3 = oo corresponds to the case Np = 2 dis- 



cussed in Sect. 2 where we found strong evidence 
for second order transition at m u d = 0. For 
mud = wis {Nf = 3), the transition is of first 
order in the chiral limit. Therefore, on the axis 
mud = 0, we have a tricritical point to* where 
the second order transition at large m 3 turns into 
first order [4]. For m s > to*, the second order 
transition line follows the m u d = axis and, for 
Til s < m*, is suggested to deviate from the ver- 
tical axis according to m u d oc (m* — to s ) 5 / 2 [7]. 
Also predicted in Ref.[4] is the MF scaling (with 
logarithmic corrections) near the tricritical point. 

Our main goal of investigations with the s 
quark is to determine the position of the phys- 
ical point in this map. I discuss the results of the 
simulations obtained so far in the next subsection. 

3.2. Simulations 

Using staggered quarks, the Columbia group 
studied this issue five years ago on a 16 3 x 4 
lattice [6] extending early studies on smaller lat- 
tices [33,34]. For the degenerate Np = 3 case, 
m u d = m s = m q , they found a first order sig- 
nal for m q a = 0.025 at (3 = 5.132. For N F = 
2 + 1, they obtained a time history suggesting 
a crossover for m u d,a = 0.025, m s a = 0.1 at 
(3 = 5.171. Their study of hadron spectrum at 
this point on a 16 3 x 24 lattice gives mjf/m p 
smaller than the experimental value suggesting 
that this m s is smaller than its physical value. 
At the same time, their large m^jmp suggests 
that their m u d is larger than the physical value. 
This implies that the physical point is located in 
the crossover region unless the second order tran- 
sition line, which has a sharp m u d dependence 
near to* (cf. Fig. 7), crosses between the physical 
point and the simulation point. 

Recently, the Tsukuba group studied this issue 
with Wilson quarks on 8 2 X 10 X 4 and 12 3 X 4 
lattices [35]. For Np = 3, they performed simula- 
tions along the K t line shown in Fig. 8 and found 
that the first order signal observed in the chiral 
limit [23] persists for (3 < 4.7 (cf. Fig. 9), while 
no clear two state signals are observed at (3 = 5.0 
and 5.5. Using a" 1 ~ 0.8 GeV for /3 < 4.7 and 
a -1 ~ 1.0 GeV for (3 = 5.0, as obtained from the 
rho meson mass extrapolated to K c , these two 
cases correspond to m q < 140 MeV and m q > 250 
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Figure 8. Phase diagram for Np = 3 with Wil- 
son quarks using the standard action [35]. For 
comparison, the K t (N t = 4) line for Np = 2 is 
included. 



MeV, respectively. (The physical s quark mass, 
giving m$ = 1.02 GeV, was shown to be about 
150 MeV with their normalization of m q .) For 
Np = 2+1, first order signals are observed for 
m ul i ~ (K u d = K c ) at both m s ~ 150 and 
400 MeV (corresponding values for K a are deter- 
mined using the a -1 values mentioned above and 
the m q a data for Np = 2-3). A recent study 
of hadron spectroscopy by the Tsukuba group for 
N F = 2 + 1 QCD at /3 = 3.5 on an 8 3 x 10 lat- 
tice shows that ~ 1.03(5) GeV at the simu- 
lation point with m s ~ 150 MeV, verifying that 
this simulation point is very close to the physical 
point in this sense. Their results on the nature 
of the transition are summarized in Fig. 10. The 
physical point clearly falls in the first order re- 
gion. 

Although both staggered and Wilson simula- 
tions give a phase structure qualitatively consis- 
tent with Fig. 7, Wilson quarks tend to give larger 
values for critical quark masses (measured by 
m^/mp etc.) than those with staggered quarks. 
This leads to the difference in the conclusions 
about the location of the physical point. How- 
ever, because both of these studies discuss that 
the deviation from the continuum limit is large at 




Figure 9. Time history of the plaquette for Np = 
3 at (3 = 4.7 on a 12 3 x 4 lattice with Wilson 
quarks using the standard action [35]. 



N t = 4, we should certainly make a calculation at 
larger N t [36] or with an improved action in order 
to draw a definite conclusion about the nature of 
the QCD transition in the real world. At present, 
increasing N t is quite painful especially with Wil- 
son quarks because of the reason discussed in 
Sect. 2.2. Therefore, the Tsukuba group applied 
their improved action to study these issues. Some 
preliminary results were reported at the confer- 
ence [29], so far with no conclusion about values 
for critical quark masses. 

4. QCD THERMODYNAMICS 

4.1. Transition temperature 

The value of the transition temperature T c in 
physical units is of great importance for phe- 
nomenological applications and for scaling tests. 

For the case of the SU(3) gauge theory, Boyd 
et al. presented new data for (3 C obtained on 
32 3 x 8 and 32 3 x 12 lattices: /3 C = 6.0609(9) 
and 6.3331(13), respectively [37]. These values 
are significantly larger than those obtained previ- 
ously on spatially small (N s /N t < 2) lattices [38] 
and remove most of the claimed discrepancy [39] 
between T c j^/a for Nt > 8 on spatially small 
lattices and that for N t = 4 and 6 using data 
for (3 C obtained on big spatial lattices [40,41] (cf. 
Fig. 11). Now we can extrapolate toward the con- 
tinuum limit much more safely: Using the infinite 
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Np= 2 pure gauge 




Figure 10. Same as Fig. 7 with the results of sim- 
ulations with Wilson quarks using the standard 
action [35]. First order signals are observed at the 
points marked by filled circles, while no clear two 
state signals are found at the points represented 
by the open circles. The values of quark mass 
in physical units are computed using a -1 ~ 0.8 
GeV for (3 < 4.7 and a" 1 ~ 1.0(1.8) GeV for 
(3 = 5.0(5.5) determined by m p (T = 0, K c ) = 770 
MeV. The real world corresponds to the point 
marked by the star. 

volume f3 c values for N t = 4 - 12 [40,41,37], Boyd 
et al. estimated 

T c Js/cj = 0.629(3), (7) 

for Nt = oo, which is much larger than the pre- 
vious estimate 0.56(3) using data for 7V t > 8 on 
small spatial lattices [39]. With sfa = 420 MeV, 
this gives T c ~ 260 MeV. 

For the case with dynamical quarks, the spatial 
lattice size is still small (N s /N t ~ 2) for N t > 6. 
This may cause a slight underestimation of T c . 
The previous status for the transition tempera- 
ture is summarized in Ref.[9]. At the conference, 
Blum reported the results of new simulations with 
Np = 2 staggered quarks performed by the MILC 
collaboration on 12 3 x 6 and 12 4 lattices [42]: 
They found that T c ~ 140 - 160 MeV^ in the 
studied range of m q a (m^/m 2 p ~ 0.1 - 0.2) both 
for Nt = 4 and 6 as can be inferred from the en- 
ergy density and the pressure plotted in Fig. 12 
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Figure 11. Transition temperature in the SU(3) 
gauge theory in units of the square root of the 
string tension. Filled circles are for the infinite 
spatial volume limit using data for (3 C obtained 
on big spatial lattices [40,41,37], and triangles are 
for previous estimates using data on small lattices 
compiled in Ref.[39]. The points at l/iV t 2 = are 
the results of extrapolations [37,39]. 

(the results for the equation of state will be dis- 
cussed in the next subsection). Here, "MeV p " 
[9] indicates that a -1 is determined by setting 
m fJ = 770 MeV at T = and at each simula- 
tion point, i.e. not in the chiral limit [43]. DeTar 
reported preliminary results of the MILC collabo- 
ration from a 24 3 X 12 lattice [44]: T c ~ 143 - 154 
MeV p for m q a = 0.008 {m\jm 2 p = 0.2 - 0.3) [44]. 
These results are consistent with the previous ob- 
servation that T c in MeV p for Np = 2 staggered 
quarks is stable, to the accuracy of about 20%, 
for a wide range of m\jra 2 p and 7V t [9]. For QCD 
with Wilson quarks, 7V t > 6 is required to get a 
stable T c [9]. No new results for iV t > 6 were 
reported this year. 

4.2. Equation of state 

Energy density e and pressure p are defined by 

p = T(j^\nZ) T (9) 

where Z is the finite temperature partition func- 
tion of the theory and V is the spatial volume in 
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T(GeV) 

Figure 12. Comparison of the equation of state 
for N F = 2 at N t = 4 (solid lines) and 6 
(dashed lines) [42]. Upper lines are for e/T 4 and 
lower lines for J>p/T 4 . The results shown are 
for m q a = 0.0125 (diamonds), 0.025 (octagons), 
and 0.1 (squares). Bursts are extrapolations to 
m q = 0. The horizontal lines give the Stephan- 
Boltzmann law for N t = 4, 6, and the contin- 
uum (lowest line). Scale for temperature T was 
fixed by m p (T = 0) = 770 MeV at each simulation 
point. 



physical units. In numerical studies on the lat- 
tice, it is useful to consider the combination 

T 9 1 7 



(10) 



where g^s are the coupling parameters of the the- 
ory [45.41]. In Eq.(10), d\nZ/dgi can be com- 
puted numerically as an expectation value of op- 
erators on the lattice, e — 3p is called "the in- 
teraction measure" because it vanishes in free 
Stefan-Boltzmann (SB) gases (see, however, the 
discussion below). On the other hand, p can be 
evaluated on the lattice using the following in- 
tegral representation derived from a thermody- 
namic relation, p = — /, with / = — TlnZ/V = 
— T91n Z/dV the free energy density for homoge- 
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Figure 13. Extrapolation to the continuum limit 
for the energy density, the entropy density and 
the pressure using data for the SU(3) gauge the- 
ory at Nt = 6 and 8 [37]. The dashed horizontal 
line shows the ideal gas limit. Separation of the 
data at T c to high and low temperature phases 
was not attempted. The hatched vertical band 
indicates the size of the discontinuity in e/T 4 [41]. 



neous systems [45]: 



pa 



3(/« 4 ) 



f9i 

-fa 4 =p(g° i )a 4 - / dg.„ 

J g « dg 



(11) 



Usually, the T = contribution is subtracted and 
gl is chosen so that the system is in the low tem- 
perature phase, to obtain p(gf) ~ 0. A conve- 
nient choice for the coupling parameter gi is (3, 
while with dynamical quarks m q a is reported to 
be useful [43]. 

For the SU(3) gauge theory, Boyd et al. studied 
e and p on 7V t = 4 - 8 lattices with large spatial 
volumes (N a /N t > 4) [37]. For the nonperturba- 
tive beta- function required in Eq.(10) and also in 
fixing the temperature scale, they used data from 
a MCRG study by the QCDTARO collaboration 
[46] and the data for f3 c discussed in Sect. 4.1. 
They extrapolated the results for p etc. to the 
continuum limit N t = oo using the volume de- 
pendence of a SB gas to leading order in 1/Nf. 
Their results show that the data for N t = 6 and 8 
are already close to their Nt = oo limit, while the 
discrepancy is significant for 7V t = 4 data (about 
20% larger in p/T 4 at T ~ 2T C ). Nevertheless, 
the previously reported sizable deviation of p/T 4 
from the SB limit for 7V t = 4 remains also in the 
N t = oo limit at T<4T C (cf. Fig. 13). 
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For the case of Np = 2 staggered quarks, Blum 
reported results for e and p by the MILC collab- 
oration on a 12 3 x 6 lattice [42] using the nonper- 
turbative beta-functions for (3 and m q a studied in 
Ref.[43]. At (3 = 5.45 and 5.53, they performed 
an extrapolation of e and p to m q = 0. Their re- 
sults, together with those for N t = 4, are shown 
in Fig. 12. Similar to the case of the pure gauge 
theory, we see that the p and e values for N t = 6 
are by 15 - 20% smaller than those for N t = 4 at 
T ~ 1.5T C . The approach to the SB limit is slow 
on both lattices. 

This slow development of p/T 4 above T c as 
well as the broad tail of the interaction measure 
at high temperatures has been considered as a 
sign for large nonperturbative interactions in the 
quark gluon plasma. 

Quite recently, Asakawa and Hatsuda [47] 
pointed out that the major part of these T depen- 
dences in the high temperature phase are caused 
by the conventional normalization p(law T) = 0, 
and the slow approach of p(T)/T 4 to p S n(T)/T 4 
by itself does not necessarily imply large non- 
perturbative effects at those temperatures: Note 
that we have a freedom to shift e and p by con- 
stants even in the case of the SB gas. These 
constant shifts are inevitable because we require 
the continuity of p at the transition between high 
and low temperature phases, as implemented in 
Eq.(ll). These shifts cause <9(l/T 4 ) deviations 
in e/T 4 and p/T 4 from the naive SB limits com- 
puted without constant terms. Because the en- 
tropy density s is determined solely from effec- 
tive degrees of freedom of the system irrespective 
of the vacuum structure etc., we do not have the 
freedom to shift s. This imposes that the shifts in 
e and p should cancel each other such that e + p 
is left untouched. Hence, the relation e = ip for 
a SB gas is also violated through this shift, and 
tails of e/T 4 , p/T 4 and (e - 3p)/T 4 in the high 
temperature phase do not reflect the width of the 
transition region directly. The width of the tran- 
sition region will be tested through the behav- 
ior of quantities defined through the derivative of 
s(T), such as the sound velocity [48] and the heat 
capacity Cy = Tds/dT. 

Asakawa and Hatsuda also showed for massless 
SB gases both in high and low temperature phases 



that e/T 4 always overshoots the naive SB value 
near T c . The behavior of e/T 4 shown in Fig. 13 
for the SU(3) gauge theory therefore suggests that 
gluons acquire an effective mass [49]. 

4.3. Gluon propagator at high tempera- 
tures 

Finite gluon mass at T > is expected also 
from a perturbative calculation. The one-loop 
gluon self-energy for SU(iVc) QCD produces a 
mass for the electric channel of the gluon propa- 
gator (the electric mass) 

m el {T) = ^(2N C + N F )/6 g(T) T (12) 

which causes the Debye screening of color electric 
charges. It is expected that the magnetic chan- 
nel of the gluon propagator also acquires a finite 
mass (the magnetic mass), although the one-loop 
perturbation theory fails to give a non-zero re- 
sult: Collecting the leading IR singularities [50], 
or by a study of a high temperature effective the- 
ory [51], we expect 

m mag (T) vO(g 2 (T)T). (13) 

These masses are important to control the in- 
frared behavior in calculations of some thermo- 
dynamic properties of the quark gluon plasma, 
such as transport coefficients. Because perturba- 
tion theory suffers from IR divergences at higher 
orders [52], a nonperturbative study is required 
to confirm these predictions. 

In the low temperature phase, however, the 
confinement makes the shape of the gluon prop- 
agator non-trivial. Previous numerical studies of 
the SU(3) gauge theory in the Landau gauge [53- 
56] reported that effective gluon masses increase 
with time separation. This implies the appear- 
ance of negative norm states. Several fits were 
attempted, so far without definite conclusions. 

This year, two studies in the high temperature 
phase were presented. Nakamura et al. simulated 
the SU(3) gauge theory at (3 = 6.8 on 48 3 x 16 and 
48 3 X 64 lattices (T/T c ~ 1.43 and 0.48, respec- 
tively) using the stochastic gauge fixing [57], and 
Heller, Karsch and Rank studied the SU(2) gauge 
theory in the Landau gauge mainly on a 32 3 x 8 
lattice at (3 = 2.6 - 3.47 (T/T c ~ 1.3 - 16) [58]. 
Both observed the increase of effective screening 
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masses of gluons with spatial separation z also in 
the high temperature phase. Nakamura et al. re- 
ported that the behavior of a gluon propagator 
is quite different at large z between the low and 
the high temperature phases, presumably due to 
the lack of confinement in the latter phase. Heller 
et al. studied the electric and magnetic channels 
separately and computed m e i and m mag in the 
high temperature phase. 2 For the magnetic mass, 
they found a good agreement with the expected T 
dependence (13), m mag (T) = 0.466(15) g 2 (T) T, 
where the two-loop formula for g(T) is used with 
A = 0.262(18) T c (g 2 (T) ~ 4.1 - 2.2 for T/T c = 2 
- 16). On the other hand, the electric mass, 
which is believed to be more reliably calculated 
[59], does not show the expected behavior (12) 
at all, even at T ~ 16T C . More work is required 
to clarify this puzzling phenomenon. Especially, 
a better understanding of the z dependence of 
gluon propagators which includes a clearer iden- 
tification of plateaus will be crucial. 

4.4. Hadronic matter below T c 

From a quenched QCD simulation with stag- 
gered quarks on a 32 3 X 8 lattice, Boyd et al. 
reported that physical quantities such as (^ , ^ r ), 
m^, f-w, info, tn p and probably tojv, do not show 
any significant T dependence up to T = 0.92T C 
[60]. Because the allowed T dependence from 
their data is much smaller than that expected pre- 
viously from calculations using the chiral pertur- 
bation theory [61] and the Nambu-Jona-Lasinio 
model [62], this makes it more difficult to ex- 
tract a phenomenologically clear evidence of a hot 
hadronic matter in heavy ion collisions. However, 
at the conference, Hioki discussed [63] that, even 
though pions etc. do not show a significant T de- 
pendence, the topological charge may have an ob- 
servable T dependence, which should affect the rj' 
meson mass [64,65]. By quenched simulations at 
(3 = 5.89 on 16 3 x 6, 16 3 x 8 and 16 3 x 16 lat- 
tices (T/T c ~ 0.93, 0.75 and 0, respectively), he 
found evidence for a T dependence of the topo- 

2 Because the effective masses are slowly increasing with 
they performed the fits by selecting data at large separa- 
tions (zT > 1) assuming that plateaus are approximately 
reached there. Propagators as a function of the momen- 
tum squared A: 2 , which were reported to be statistically 
stable [55,56], were not studied. 



logical charge distribution. This information can 
be used to compute m v i. His results suggest that 
there exists an observable T dependence in m v i. 

4.5. Valence quark chiral condensate 

The Columbia group extended their previous 
study [66] of valence quark mass dependence in 
Np = 2 QCD with staggered quarks using 16 3 x 4 
and 32 3 x 4 lattices [67]. They studied the va- 
lence quark chiral condensate (CO as a function 
of the valence quark mass m va i and the sea quark 
mass m q . (££) essentially measures the density 
of Dirac eigenstates with eigenvalues of 0(m va \). 
It can be shown that {CC)( m vai = 0,m 9 ) = is 
a sufficient condition for the restoration of chi- 
ral symmetry. Their data show that m va i depen- 
dence of (££) is actually sensitive to the phase. 
They also presented results for ((J dx-ip^Ti-ip) 2 ) — 
(( J dxtjiTiij}) 2 ) in the vicinity of T c which sug- 
gest that the axial Ua(1) symmetry is broken just 
above T c . 

5. FINITE CHEMICAL POTENTIAL 

In spite of phenomenological urgency, QCD at 
finite chemical potential fi on the lattice [68] re- 
mains to be a challenging issue because of the 
complex action problem [8]. The standard argu- 
ment suggests that, at T = 0, the chiral symme- 
try is restored at fi c « tojv/3 where mjv is the 
nucleon mass [69]. 

On the other hand, quenched QCD simula- 
tions, which ignore the fermion determinant con- 
taining the complex phase, show pathological be- 
havior at m v /2<[i<m N /3 for T = [69-71]: 
Finite baryon number and energy densities ap- 
pear at /i>m T /2, and ( 1 l r \E r ), which is expected 
to be independent of fi below fi c , begins to de- 
crease at fi ~ m T /2. For /i>mjv/3, a re- 
cent study at (3 = shows that all observables 
have their limiting vales: ($\f r ) = etc. [71]. 
Simulations for m T /2 < fi < tojv/3 require a large 
amount of the CPU time due to large fluctua- 
tions in observables and slow convergence of the 
algorithm used to invert fermion matrices. To un- 
derstand this unexpected behavior in the interval 
<[i< tojv /3, the effect of the quenching ap- 
proximation [72], artifacts from finite lattices [73], 
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and problems due to the use of staggered quarks 
[74] are examined, so far, however, without a def- 
inite conclusion. 

In order to see if this unexpected behavior is 
caused by a physical transition at fi « m T /2, it is 
important to clarify the role of dynamical pions at 
finite chemical potential [70,71]. Recently, Hands 
et al. performed a full simulation of the 2+1 D 
Gross-Neveu model with staggered fermions at 
fi > [75]. Although the model lacks the fea- 
ture of confinement, it contains a massless pion 
at fj, < [i c in the chiral limit as a mesonic bound 
state of elementary fermions. An important point 
is that the fermion determinant of the theory is 
real and positive semidefinite also at fi > 0. The 
results of their simulation are completely consis- 
tent with naive expectations and a 1/Np calcu- 
lation: no unexpected behavior is observed at 
fi ~ m T /2. This suggests that the pathological 
behavior of quenched QCD at < fi < tojv/3 

is not physical. Although other possibilities are 
not excluded [73,74], it is natural to expect that 
a full QCD simulation will lead us to a simple 
phase structure. At present, full QCD simula- 
tions are attempted only on small lattices (typi- 
cally 4 4 ) and/or in the strong coupling limit [76- 
80]. More efforts are needed to solve the problem 
of complex actions. 3 

6. CONCLUSIONS AND PROSPECTS 

I reviewed recent progress in finite temperature 
lattice QCD. We found evidences that the scal- 
ing properties near the chiral transition in two 
flavor QCD are consistent with the 0(4) scaling 
both with staggered and Wilson quarks. It was 
pointed out that use of some improved action is 
essential to obtain the scaling behavior for Wil- 
son quarks on small N t lattices (presumably for 
N t < 18), by eliminating lattice artifacts with the 
standard action. Also a study on the structure of 
the critical line at T > with Wilson quarks was 
reported. These progresses provide us with more 
rigid conceptual and technical bases for Wilson 
quarks. For more confirmation of the 0(4) scal- 
ing, further studies and better accuracy of data 

3 A different approach, the static approximation, was pre- 
sented at the conference [81]. 



are required: For staggered quarks, data on larger 
lattices and at smaller quark masses are needed 
to check consistency with the 0(4) scaling, and, 
in future, to distinguish between 0(4) and 0(2) 
exponents. For Wilson quarks, a direct determi- 
nation of critical exponents should be attempted. 

I also reviewed the status of investigations in 
QCD including the strange quark. Although both 
staggered and Wilson simulations performed on 
N t = 4 lattices give a phase structure qualita- 
tively consistent with the expected one shown in 
Fig. 7, Wilson quarks tend to give larger criti- 
cal quark masses in comparison with staggered 
quarks. This consequently leads to the difference 
in the conclusion about the location of the real 
world in Fig. 7. Because the deviation from the 
continuum limit is large in the both simulations 
at 7V t = 4, we should certainly make a calculation 
at larger N t or with an improved action in order 
to draw a definite conclusion about the nature of 
the transition in the real world. 

I finally reviewed the progress in other topics 
of finite temperature/density QCD. The results 
are both encouraging and challenging and it is 
important to perform further studies in these di- 
rections. Improvement of the lattice action may 
prove to be useful also in these studies. 

I am indebted to S. Aoki, T. Blum, S. Chan- 
drasekharan, N. Christ, T. Hatsuda, S. Hioki, 
F. Karsch, A. Kocic, J.B. Kogut, E. Laermann 
A. Nakamura, D.K. Sinclair, and A. Ukawa for 
valuable discussions and providing me with orig- 
inal data. I am also grateful to my collabora- 
tors, Y. Iwasaki, S. Kaya, S. Sakai, and T. Yoshie 
for their support and encouragements. Finally I 
would like to thank Y. Iwasaki, A. Ukawa and W. 
Bock for critical comments and suggestions on the 
manuscript. This work is in part supported by 
the Grant-in- Aid of Ministry of Education, Sci- 
ence and Culture (Nos.07NP0401 and 07640376). 
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